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Abstract. In this letter we present the results of a wavelet analysis of the radio light curve of the BLLac Object 
PKS B0048— 097 from the University of Michigan Radio Astronomy Observatory monitoring program at 8 GHz 
during twenty-five years, from 1979 to 2004. The results show a remarkable periodicity of 450-470 days in the 
early 1980s that changed to a ~585 day periodicity in the late 1980s to early 1990s. A less pronounced ~400 day 
periodicity is found after ~ 1995. Very-long-baseline interferometry imaging at 15 GHz shows dramatic structural 
changes in the usually unresolved source between two epochs, 1995.57 and 2002.38. The pronounced northward 
directed jet seen in the 2002 image differs by more than 90° in direction from the source structure found in earlier 
epochs. These findings make PKS B0048— 097 a primary target for multi-wavelength observations and intensive 
radio monitoring to decipher the blazar-variability phenomenon. 
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1. Introduction 

Combined multi- frequency (radio to X-ray/7-ray) obser- 
vations of active galactic nuclei (AGN) can yield impor- 
tant insights into the physics of relativistic outflows as- 
sociated with super-massive black holes. In particular the 
emission of blazars (BL Lac objects and OVV quasars, see 
Ulrich, Maraschi & Urry 119971 for a description on the 
nature of blazars) is believed to be dominated over the 
whole electromagnetic spectrum by the most compact re- 
gions of relativistic jets. Blazars exhibiting quasi-periodic 
behaviour in principle allow the inter-relation between the 
compact radio jet and the broadband spectral energy dis- 
tribution to be deciphered via quasi-simultaneous multi- 
frequency observations, especially if the time scales are 
small enough to consider variability data over multiple 
cycles, but such objects are rare. 

Hughes, Aller, & Aller l|1998|l find a persistent mod- 
ulation of the total flux and polarisation of the BLLac 
object OJ287. The dominating period of ~1.12yr in the 
1980s was later displaced by a strong ~1.66yr periodicity 
during the 1990s. The relationship between these two vari- 
ations was interpreted by Hughes, Aller, & Aller 



optical variability of OJ 287 (period between major out- 
bursts ~ 11.6 yr), has been interpreted, e.g., in terms of a 
"precessing-jet" model by Abraham l|2UUU|l . Stirling et al. 
( 2003) suggest an oscillating "nozzle" structure of the in- 
ner jet of BL Lac with a period of ~ 2 years, based on radio 
millimetre flux-density monitoring and VLBI observations 
between 1998 and 2001 (see also Mutel & Denn EES). 
Villata et al. (2004) analyse the optical and radio long- 
term variability of BL Lac and find a gradually lengthen- 
ing ~ 8 yr periodicity. The optical and radio long-term 
variability of the BLLac object AO 0235+16 has been in- 
vestigated by Raiteri et al. (|2001[1 revealing a possible 5-6 
year periodicity based on long-term monitoring data, in 
the radio regime particularly on data of the UMRAO 1 
database. Ostorero, Villata & Raiteri ( 2004) applied the 
helical jet model of Villata & Raiteri i 1 !!!)!! i to the case 
of AO 0235+16, interpreting the data in terms of Doppler 
factor variations due to changes of the viewing angle be- 
cause of the helix rotation. In all these cases, attempts 
have been made to understand the periodic behaviour as 
a result of orbital motion of two black holes in a binary 



terms of a "shock-in-jet" model while the complex periodic 
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system, helical jet structures, shocks, and instabilities of 
the disk or jet-plasma flow. 

In this letter, we report on the detection of a par- 
ticularly strong, and persistent quasi-periodic long-term 
modulation of the compact-jet emission of the BL Lac ob- 
ject PKS B0048— 097. The comparatively short time scale 
of ^350 days to ^600 days and the large amplitude of 
the variability make PKSB0048-097 a particularly well- 
suited source for future multi-frequency campaigns on the 
broadband spectral emission of blazar sources and tests of 
blazar-periodicity models. 

In Sect. 2, we present the observational data and their 
analysis. We discuss the results of a cross- wavelet anal- 
ysis of twenty five years of the University of Michigan 
Radio Astronomy Observatory (UMRAO) monitoring of 
PKSB0048— 097 and discuss its parsec-scale structure 
during two epochs observed as part of the Very Long 
Baseline Array (VLBA) 2 cm Survey. Finally, we summa- 
rize our results and their implications for future blazar 
studies, in Sect. 3. 

2. Background, Observations and Analysis 

PKSB0048-097 is a BLLac object with unknown red- 
shift. A lower limit of z = 0.5 is postulated by 
Falomo ifTOMjl from HST data. The flux density of 
PKS B0048-097 has been measured since 1979 at 4.8 GHz, 
8.0 GHz, and 14.5 GHz using the University of Michigan 
26 m paraboloid. The source is observed once every three 
months at all three frequencies as part of the UMRAO 
BLLac Observing Program (Aller et al. 119991 Aller, Aller 
& Hughes 2003- Its light curve is best sampled at 8.0 GHz 
(see Fig. ^| exhibiting pronounced maxima during several 
well defined time periods. 

2.1. Wavelet Analysis 

We used a continuous wavelet analysis, which quantifies 
the behaviour of a signal on different temporal scales, to 
investigate the available UMRAO data for periodicities. 
As a function of time, the signal is convolved with a lo- 
calised wave-packet, that is translated along the series, 
for a number of 'dilations' of the wave-packet (Hughes et 
al. I1998|l . This technique has the great advantage of pre- 
serving temporal locality: a gap in the time series will be 
evident along the corresponding line in transform space, 
and events that are distinct in the signal will have dis- 
tinct counterparts in transform space. A Morlet wavelet 
is particularly suited to the analysis of time series, as it 
is complex, so the real part of the transform exhibits an 
oscillatory behaviour corresponding to periodicity in the 
time series, while the modulus provides a measure of the 
power in different components of the signal. 

Such an analysis has been applied to the time series 
for PKS B0048— 097, and, as seen in Fig.^ a distinct pat- 
tern is seen in both the real part of the source wavelet 
transform, and its modulus, corresponding to a quasi- 
periodic component containing a modest fraction of the 
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Fig. 1. A continuous Morlet wavelet transform of the time 
series data of the flux density for PKS B0048-097 taken at 
the UMRAO at 8.0 GHz. The panels show the flux density 
time series (top), the real part of the transform (middle) 
and the modulus (power) of the transform (bottom). The 
real part displays clearly the signature of a quasi-periodic 
component with a slowly varying ~ 1-2 yr period. Power 
at this dilation extends across the entire series. The hor- 
izontal lines in the middle and bottom panel enclose the 
95 % significance region of the detected periodicity deter- 
mined by a cross- wavelet transform analysis (see text). 
The hatched areas of the transform panels show the cone 
of influence, where edge effects resulting from the finite 
length of the time scries decrease the wavelet coefficient 
power by more than a factor e 2 . 



overall power, and masked in a Fourier power spectrum, 
because of drifts in time scale across the data window. In 
the cross wavelet transform technique, described, e.g., by 
Kelly et al. ( 2003 ) , the continuous transform of the signal 
is convolved with a set of template periodic signals. This 
analysis has been used to quantify the result (apparent on 
visual inspection of the top panel at Fig.^l, establishing 
a time scale of 385-470 days with a confidence of > 99.9% 
that this component does not arise by chance from random 
patterns associated with a lag-1 autoregressive process. 

The varying time scale of the disclosed periodicity 
has been further investigated using a Lomb-Scargle pe- 
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Fig. 2. Lomb-Scargle periodogram for three different 
parts of the time series of UMRAO PKS B0048-097 8 GHz 
monitoring data. The dotted and dashed vertical lines 
show the 95% and 99.9% significance lines, respectively, 
derived from the cross wavelet transform of the full light 
curve. 



riodogram analysis (Tomb IT9761 Scargle l*L982[l . In a L-S 
periodogram, the classical discrete Fourier periodogram is 
redefined in such a manner to as to make it invariant to a 
shift of the origin of time. The L-S periodogram analysis 
has been performed using the PERIOD software package 2 . 
The L-S analysis between November 1979 and September 
1985 shows maximal power at ~ 400-460 days and a sec- 
ondary signal at ~ 580 days (see Fig.0). Between October 
1985 and February 1995, only one narrow peak at ~ 585 
days is present. After October 1995, the highest power is 
found at a time scale of ~ 405 days, that is considerably 
less pronounced than the periodicities at the earlier time 
ranges. The significance of the individual peaks in Fig.[5]is 
difficult to quantify because of the not a-priori known fre- 
quency dependence of the underlying noise process 3 . We 
base our detection on the highly-significant signal in the 
global power spectrum. Figure El illustrates that over the 
full ~ 25 years of UMRAO 8 GHz radio light curve mon- 
itoring of PKSB0048— 097 the maximal power found by 
the L-S analysis is located within the significance ranges 
found by the cross- wavelet analysis, but that its time scale 
and power vary with time. 



2 The period software package is provided by the Starlink 
Project which was run by CCLRC on behalf of PPARC. See 
: / / star-www .rl.ac.uk 
a Only the assumption of a white noise process would al- 
low a simple significance threshold parallel to the frequency 
axis to be placed. Admitting a red noise process would require 
us to consider a range of possible power laws for the process, 
which would lead to a range of significance areas (not parallel 
to the frequency axis). Note, that the formal L-S false-alarm- 
probability lies below 1 % (with 95 % confidence) for all major 
peaks in Fig. 



2.2. Structural Variability 

On parsec scales, Shen et al. I|1997fl and Gabuzda, 
Pushkarev & Cawthorne ( 199j||) both report a core-jet 
morphology of PKS B0048— 097, however the reported jets 
differ by ^40 degrees in position angle (P.A.). While Shen 
et al. I|1997|) report a P. A. of ~ -160° (south-westward) in 
epoch 1992.9 from 5 GHz VLBI observations, Gabuzda et 
al. lfTU9"9"|l find a jet at P.A. ~ +160° (south-eastward) 
in epoch 1992.2 also from 5 GHz data. In both cases 
the (it, v)-coverage was not optimum. In the VLBA 2 cm 
Survey observations before 2002, PKSB0048-097 showed 
no clear resolved structure. FigureOlshows the milliarcsec- 
ond structure of PKSB0048-097 in epochs 1995.57 and 
2002.38 (compare Table [TJ. In 1995.57, a weak westward 
directed jet is found, partially resolved only by the longest 
east-west baselines. The (u,v)-d&t& obtained in 2002.38 
show clear evidence for a resolved north-south structure 
along P.A. ~ —30° with the jet pointing in a direction 
more than 90° different from what was reported by Shen 
et al. 119971 and Gabuzda et al. l|1999fl . Considering these 
discrepancies between different works, PKSB0048— 097 
might represent a case of extreme jet-ejection-angle vari- 
ation. This would be of particular interest in conjunction 
with a putative periodicity of the radio light curve as re- 
vealed from the UMRAO data above. In the scenario of a 
precessing jet, PKS B0048— 097 might represent a highly 
attractive target to study the broadband jet emission of a 
BLLac object at different angles to the line of sight. 

On the basis of the available data, it is not possible to 
derive the time scales of the structural variability. VLBA 
2 cm Survey observations of PKS B0048— 097 have been 
conducted during two additional epochs between 1995 and 
2002 but the data quality in both cases does not allow us to 
compare the milliarcsecond structure directly to the two 
images presented here, due to shorter integration times 
and sub-optimal (u, w)-coverage. After 2002, the source 
has been observed several times as part of the VLBA 
2 cm Survey continuation project, MOJAVE. An analysis 
of these observations, which provide full polarimetric in- 
formation, is currently being performed. A first MOJAVE 
image of PKSB0048-097 in 2003 showing a northward- 
directed jet similar to the 2002 source structure can be 
found in Lister & Homan ((,2005 ). 

3. Summary and Implications 

The BLLac object PKSB0048-097 exhibits strong struc- 
tural variability on sub-milliarcsecond scales and pro- 
nounced radio, optical (see, e.g., Pica et al. 119881) and 
X-ray flux variability. We have demonstrated that the ra- 
dio light curve of PKSB0048— 097 has a strong modula- 
tion with a period of ~ 350 days to ~ 600 days at 8 GHz. 
Additional insight may be gained from future polarimet- 
ric studies which have not been considered in this letter. 
From all wavelet transforms of sources monitored as part 
of the UMRAO program, only for OJ 287 a comparably 
strong persistent signal has been found. The strength of 
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Fig. 3. VLBI results of PKS B0048-097 at 15 GHz. Top 
Panels: Naturally weighted images of PKSB0048-097 in 
1995.58 (left) and in 2002.38 (right) at A2cm obtained 
from VLB A 2 cm Survey observations. Image parameters 
are given in Table A jet towards the north is visible 
in the 2002 image while a slight extention to the west 
is seen in the earlier 1995 image. Best fitting Gaussian- 
modcl-component approximations to the (u, u)-data are 
indicated as ellipses. Middle and bottom panels: The pro- 
jected visibility along P.A. 160° and —70° for both epochs. 
Clear deviations from a flat point-source visibility indicate 
an elongation of the source along different position angles. 

Table 1. VLBI Results 





Image parameters a 


Model fit parameters 11 


Epoch 


Beam L Speak rms e 


Comp. 


r 8 P.A." 


S' 




[(masxmas), ] [Jy/beam] [mjy/beam] 




[mas] [°] 


[Jy/beam] 


1995-07-28 


(1.3 x 0.5); -5.2 1.21 0.3 


Core 





1.186 






Wl 


0.2 -97.4 


0.052 


2002-05-18 


(1.3 x 0.5); -2.6 1.09 0.2 


Core 





1.074 






Nl 


0.6 -22.2 


0.110 






N2 


2.1 8.7 


0.012 



a The lowest contour is 1 mjy/beam in both images, contours increase 
by factors of v^2; b Formal uncertainties in the model fits are too 
small and therefore are not printed in the table; c Size and orientation 
of the restoring beam; d Peak flux density per beam area; c Noise 
level; f Individual Gaussian model-fit component; s Distance from core 
component; h Position angle of model-fit component; 1 Flux density of 
model-fit component per beam area. 



the modulation and the comparatively short time scale 
makes PKSB0048— 097 a prime object for coordinated 
broadband observing campaigns. In particular, the corre- 
lation of the varying VLBI-jet ejection angle of the source 
with radio (and higher energy) light-curve evolution pro- 



vides a tool to test precessing-jet models against alter- 
native models (e.g., intrinsic instabilities of the flow) of 
periodic light-curve modulation in blazars. Unlike most 
other sources, the short 350 day to 600 day time scale of 
PKS B0048-097 allows these alternatives to be tested ob- 
servationally in a relatively short time. 
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